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Abstract 

The momentum transfer dependence of the LSP-nucleus elastic scattering cross sec¬ 
tions is studied. New imput SUSY parameters obtained in a phenomenologically allowed 
parameter space are used to calculate the coherent rate for various nuclear systems and 
the spin matrix elements for the proposed Pb target. The results are compared to those 
obtained from other cold dark matter detection targets. 
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1 Introduction 


Recently, the possibility to directly detect the lightest supersymmetric particle (LSP), a can¬ 
didate for cold dark matter, via the recoiling of a nucleus in the elastic scattering process: 

X + (A, Z) ^ X + (A,Zr ( 1 ) 

(x denotes the LSP) has been proposed. |I|, In the present study we proceed as follows: 

1) We write down the effective Lagrangian at the elementary particle (quark) level obtained 
in the framework of supersymmetry as described in refs. |]l|, ^ 

2) We go from the quark to the nucleon level using an appropriate quark model for the 
nucleon. Special attention in this step is paid to the scalar couplings, which dominate the 
coherent part of the cross section and the isoscalar axial currents which depend on the assumed 
quark model. 

3) We compute the relevant nuclear matrix elements using as reliable as possible many boby 
nuclear wave function hopping that, by putting as accurate nuclear physics input as possible, 
we will be able to constrain the SUSY parameters as much as possible. 

4) We calculate the modulation of the cross sections due to the earth’s revolution around 
the sun and the motion of the sun around the center of the galactic disc by a folding procedure. 

Our goal is to calculate LSP-nucleus event rates for Pb nucleus by using new imput SUSY 
parameters ii obtained in a phenomenologically allowed parameter space. We focus on the 
spin matrix elements of since this target, in addition to its experimental qualihcations, 

has the advantage of a rather simple nuclear structure. Furthermore, its spin matrix element, 
especially the isoscalar one, does not exhibit large quenching as that of the light and up to now 
much studied and nuclei. We compare our results to those given from other cold 
dark matter detection targets. 

The total cross section of the LSP-nucleus reaction (1) can be written as 


a 


nip (1 -I- r]Y A 
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( 2 ) 


where mi, nip is the LSP, proton masses, rj = mi/rUpA and cxo — 0.77 x 10“^®cm^. The 
parameters /g, fy, /^, with r = 0,1 an isospin index, describe the scalar, vector and axial 
vector couplings and are determined in various SUSY models p, |] (see ref. 0]). The momentum 
transfer enters via u as 
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(b is the harmonic oscillator parameter). In Eq. (|^), the terms in the square brackets describe 
the coherent cross section and involve the integrals Ip{u) [p = 0,1) while the other terms 
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Figure 1: (a) Plot of the integrals and Io{u) for Pb, where u is given by Eq. (3). We see 

that Ill is quite a bit less retarded compared to Iq. (b) Kq integrals (for 1=0 and 1=1) entering 
the dominant scalar part of the event rate. 


describe the spin dependence of the cross section by means of the integrals Ipp'{u) (p, p' = 0,1, 
isospin indices) which are normalized so as /p(0) = 1 and /pp'(O) = 1. For their definitions see 
ref. where they are calculated by using realistic nuclear form factors. 

In Fig. 1(a) we show the variation of Iq{u) (the others behave similarly). We see that u 
can be quite big for large mass of the LSP and heavy nuclei even though the energy transfer 
is small (< lOOiFeE). The total cross section can in such instances be reduced by a factor of 
about five. 

The spin matrix element of heavy nuclei like Pb has not been previously evaluated, since 
one expects the relative importance of the spin versus the coherent mode to be more pronounced 
on light nuclei. However, the spin matrix element in the light isotopes is quenched, while that 
of does not show large quenching. For this feature, we recently proposed Pb nucleus as 
an important candidate in the LSP detection. Furthermore, this nucleus has some additional 
advantages as: i) it is believed to have simple structure (one 2pl/2 neutron hole outside the 
doubly magic nucleus ^°®P6) and ii) it has low angular momentum and therefore only two 
multipoles A = 0 and A = 2 with a J-rank of k = 1 can contribute even at large momentum 
transfers. 

The dependence of the cross section on the momentum transfer for ^°^P6 has been exten¬ 
sively studied in ref. In Fig. 1(a) the variation of the integral /ii(m), which describe the spin 
part of ^°^P6, is shown. 


2 Folding of the cross section with a Maxwellian distri¬ 
bution 


Due to the revolution of the earth around the sun, the event rate becomes modulated. This 
effect is studied by folding the total cross section Eq. (^ with an appropriate distribution. If 
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we assume a Maxwell type distribution, which is consistent with the velocity distribution of 
LSP into the galactic halo, the counting rate for a target with mass m takes the form [Q 


/ dN 
' dt 


p(0) m , -r— 

-V< > < S > 

mi Amp 


where p(0) = 0.3GeV"/cm^, the LSP density in our vicinity and < S > is given by 


( 4 ) 
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The parameters Jq, Jp, Jp^j describe the scalar, vector and spin part of the velocity averaged 
counting rate, respectively. They are functions of Pq = vq/c^ 6 = 2v\/vo and mq, where vq is 
the velocity of the sun around the galaxy, vi the velocity of the Earth around the sun and 
Mo is given by an expression like that of Eq. (||) with f3 —^ /do- Since 5 ~ 0.27 << 1, we 
can expand the J-integrals in powers of 5 and retain terms up to linear in 5. Thus, for each 
mechanism (vector, scalar, spin) we obtain two integrals associated with I = 0 and I = 1. The 
most important integrals are shown in Fig. 1(b). For some others see ref. By exploiting 
the above expansion of K-integrals, the counting rate can be written in the form 


/ dN \ / dN \ 

\~dr/ ^ \ir/o 


(1 + h cosa) 


( 6 ) 


where a is the phase of the earth’s orbital motion and {dN/dt)^ is the rate obtained from the 
I = 0 multipole, h represent the amplitude of the oscillation, i.e. the ratio of the component of 
the multipole / = 1 to that of the multipole / = 0. 


3 Results and discussion 


The three basic ingredients of our calculation were the input SUSY parameters, a quark model 
for the nucleon and the structure of the nuclei involved. The input SUSY parameters used 
have been calculated in a phenomenologically allowed parameter space (cases #1, #2, ^3) as 
explained in ref. 

In Tables 1, we compare the spin matrix elements at g = 0 for the most popular targets 
considered for LSP detection and We see that, the spin matrix elements 

of are even a factor of three smaller than those for obtained in ref. ||^ 

The spin contribution, arising from the axial current, was computed in the case of 
system (see Tables 2, 3 and 4). For the isovector axial coupling the transition from the quark to 
the nucleon level is trivial (a factor of qa = 1.25). For the isoscalar axial current we considered 
two possibilities depending on the portion of the nucleon spin which is attributed to the quarks. 
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Table 1: Comparison of the static spin matrix elements for three typical nuclei. 


Component 


''^Gegj2+ 

^“S'q/2+ 


0.231 

1.005 

0.204 

Pi(0)Po(0) 

-0.266 

-1.078 

-0.202 

P§(0) 

0.305 

1.157 

0.201 


Table 2: The spin contribution in the LSP —Pb scattering for two cases: EMC data and 
NQM Model. The LSP mass is mi = 126, 27, 102 GeV for #1, #2, #3 respectively. 



EMC DATA NQM MODEL 

Solution 

< dN/dt >0 (y-^Kg-^) h 

< dN/dt >0 {y ^Kg h 

#1 

#2 

#3 

0.285 X 10-^ 0.014 

0.041 0.046 

0.012 0.016 

0.137 X 10-^ 

0.384 X 10-2 0.056 

0.764 X 10-2 0.017 


indicated by EMC and NQM. The ground state wave function of obtained by 

diagonalizing the nuclear Hamiltonian in a 2h-lp space which is standard for this doubly magic 
nucleus. 

For the coherent part (scalar and vector) we used realistic nuclear form factors. In ref. P| 
we studied three nuclei, representaves of the light, medium and heavy nuclear isotopes (Ca, Ge 
and Ph) for three different quark models. In table 4 we present the results for Pb obtained for 
these models denoted by A (only quarks u and d in the nucleon) and B and C (heavy quarks 
included) for the LSP masses shown in table 3 (cases 7)^4-#9). We see that the results vary 
substantially and are sensitive to the presence of quarks other than u and d into the nucleon. 

The total cross section is almost the same for LSP masses around 100 GeV. We obtained 
results in the context of the quark models NQM, EMC, for SUSY models #l-#3 0 (Tables 2 
and 3) and SUSY models #4-#9 [Q (Tables 4 ). 


4 Conclusions 


In the present study we found that for heavy LSP and heavy nuclei the results are sensitive 
to the momentum transfer as well as to the LSP mass and other SUSY parameters. From the 
Tables 3 and 4 we see that, the results are also sensitive to the quark structure of the nucleon 
in the following sense: i) The coherent scalar (associated with Higgs exchange) for model A (u 
and d quarks only) is comparable to the velocity suppressed vector coherent contribution. Both 
are at present undetectable, ii) For models B and C (heavy quarks in the nucleon) the coherent 
scalar contribution is dominant. Detectable rates < dN/dt >o> 100 y~^Kg~^ are possible in a 
number of models with light LSP. 

The spin contribution is sensitive to the nuclear structure. It is undetectable if the LSP is 
primarily a gaugino. 

The folding of the event rate with the velocity distribution provides the modulation effect 
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Table 3: Ratio of spin contribution Ph/'^^Ge) at the relevant momentum transfer with the 
kinematical factor 1/(1 + rff, rj = mi/Arrip. 


Solution 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

#9 

mi (GeV) 

126 

27 

102 

80 

124 

58 

34 

35 

50 

NQM 

0.834 

0.335 

0.589 

0.394 

0.537 

0.365 

0.346 

0.337 

0.417 

EMC 

0.645 

0.345 

0.602 

0.499 

0.602 

0.263 

0.341 

0.383 

0.479 


Table 4: Event rates for Ph. The LSP mass is mi = 80, 124, 58, 34, 35, 50 GeV for the cases 
4/4: — respectively. Cases #9 are no-scale models. The values of < dN/dt >o for Model 
A and the Vector part must be multiplied by xl0“^. 



Scalar Part 

Vector 

Part 

Spin Part 



/M\ 

\ dt / 0 


h 

~TIK\ 

\ dt /i\ 

h 


tmT) 

\ d.t / 

h 


A 

B 

C 



EMC 


NQM 


#4 

0.03 

22.9 

8.5 

0.003 

0.04 

0.054 

0.8010- 


0.1610-^ 

0.015 

#5 

0.46 

1.8 

1.4 

-0.003 

0.03 

0.053 

0.3710- 

-3 

0.9110-3 

0.014 

#6 

0.16 

5.7 

4.8 

0.007 

0.11 

0.057 

0.4410- 

-3 

0.1110-2 

0.033 

#7 

4.30 

110.0 

135.0 

0.020 

0.94 

0.065 

0.67 


0.87 

0.055 

#8 

2.90 

73.1 

79.8 

0.020 

0.40 

0.065 

0.22 


0.35 

0.055 

#9 

2.90 

1.6 

1.7 

0.009 

0.95 

0.059 

0.29 


0.37 

0.035 


h. In all cases it is small, less than ±5%. 
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